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We use local and global magnetometry measurements to study the influence of magnetic domain 
width w on the domain-induced vortex pinning in superconducting/ferromagnetic bilayers, built of 
a Nb film and a ferromagnetic Co/Pt multilayer with perpendicular magnetic anisotropy, with an 
insulating layer to eliminate proximity effect. The quasi-periodic domain patterns with different 
and systematically adjustable width w, as acquired by a special demagnetization procedure, exert 
tunable vortex pinning on a superconducting layer. The largest enhancement of vortex pinning, by 
a factor of more than 10, occurs when w ~ 0.31/im is close to the magnetic penetration depth. 

PACS numbers: 74.25Ha, 74.25Qt, 74.78Db, 74.78Fk 
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One of the most important parameters for the appli- 
cation of type II superconductors is the critical current 
density, Jc- To achieve high Jc it is necessary to pin vor- 
tices, which exist in the mixed state. The usual method 
for pinning vortices utilizes microscopic defects in the ma- 
terial that locally suppress superconductivity and trap 
normal vortex cores. Magnetic pinning (MP) instead re- 
lies on the electromagnetic interaction between the vor- 
tex magnetic field and the stray field generated by the 
magnetic texture in the vicinity of the superconductor 
surface [l|-|3[. Since magnetic pinning acts on a length 
scale comparable to the penetration depth, A, it is domi- 
nant at temperatures close to the superconducting tran- 
sition temperature, T^, where pinning by defects becomes 
ineffective. 

A novel magnetic pinning method utilizes a planar 
superconductor(S)/ferromagnet(F) bilayer (SFB) sepa- 
rated by a thin insulating layer, which eliminates proxim- 
ity effects [ll-Q- ^s previously suggested, the magnetic 
domains in the F layer create pinning centers for the vor- 
tices in the S layer so that tuning the domain structure 
can result in tuning Jc \%\^. While the principle of 
magnetic pinning in the SFB's has been extensively dis- 
cussed theoretically |2h11| and demonstrated experimen- 
tally |12h26I | , the magnitude of the pinning enhancement 
thus far reported has been modest with a factor of no 
more than 3, and sometimes even suppression of pinning 
has been described [18|. Furthermore, it has not been fea- 
sible to compare various magnetic pinning results, since 
the different SFB bilayers used in these studies would 
render such comparison impractical. 

In this work we show that it is possible to systemati- 
cally vary and achieve much higher enhancement of Jc in 
planar SFB's with Nb as the S layer and Co/Pt multi- 
layer with perpendicular magnetic anisotropy (PMA) as 
the F layer. The PMA in the F layer assures a direct ef- 
fect in magnetic pinning that is not well defined in F layer 
with in-plane anisotropy. Using a single SFB, we use a 
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FIG. 1: (Color online)(a) Schematic of the samples A and 
B. (b) Hysteresis loop for sample B at T = 10 K. (c-e) MFM 
images (10 fim x 10 fim) at 300 K for sample A demagnetized 
at angles 61 = 0°, 20° and 90°. (f) w vs 9 for samples A, B 
and C (lines are guides to the eye). 



special demagnetization method to continuously tune the 
width of the domains, w, in the quasi-periodic stripe do- 
main pattern of the F layer with equal +/— domain, and 
observe the dramatic effect of tuning on vortex pinning 
using the global (SQUID) and the local (Hall sensors) 
magnetometry measurements. Using a single SFB with 
an identical S layer, we have quantitatively determined 
Jc enhancement as a result of changing domain width. 
We have observed Jc enhancement in excess of 10. 

The SFB's were grown by sputtering at room tem- 
perature on Si(lOO) substrates with the same sequence 
of Si(10)/Pt(10)/[Co(0.6)/Pt(0.3)]8/Si(10)/Nb(t), with 
the thickness denoted in nanometers in parentheses (Fig. 
1(a)). The Si(lO) layer between [Co/Pt]8 and Nb elimi- 
nates the proximity effect. Two samples (A and B) with 
nominally t — 76 nm were studied by magnetometry, 
whereas the third sample (C) with t = 20 nm with an 
additional Si(lO) protective layer on the top was used to 
extract the vortex activation energy from transport ex- 
periment, as described separately [23|. All three samples 



(A, B, C) display a square hysteresis loop in the normal 
state, with the coercive fields at 10 K oi He = 720 Oe 
(A), 775 Oe (B), and 750 Oc (C) (Fig. 1(b)). From 
magnetoresistance measurements of SFB magnetized to 
saturation we extract the superconducting parameters of 
the Nb layer of T^ = 8.56K (A), 8.42 K (B), and 7.95 K 
(C), coherence length, ^(0) = 14.5 nm (A), 12.4 nm (B), 
and 12 nm (C), and estimated Ginzburg-Landau param- 
eter for dirty-limit superconductors [27| of k « 3.8 (A), 
5.2 (B), and 5.6 (C). 

A well known method to acquire multi-domains with 
equal +/— domains in an F layer with PMA is AC de- 
magnetization with field perpendicular to the film plane. 
However, we can also obtain domain pattern with equal 
+/— domains but with different width by demagnetiza- 
tion with an AC magnetic field at an angle with the 
sample plane. Figs. 1(c)- (e) show MFM images of sam- 
ple A at 300 K for various 6''s. Indeed, the images reveal 
equal amount of positive and negative domains, but the 
average width w increases with increasing 9. From these 
MFM images, using two-dimensional Fourier analysis we 
extract the mean value of w and the standard deviation, 
as shown in Fig. 1(f). Demagnetization a.t 9 — 0° re- 
sults in a similarly small w « 0.31 - 0.35 /im for all sam- 
ples, but the width w increases with 9 with a more rapid 
growth of w at small 9. This increase of w with 9 is caused 
by the fact that in ferromagnets with PMA the mag- 
netic moments are out-of-plane in the uniform domain 
area but with an in-plane component within the domain 
walls. With increasing 9 the demagnetization field has a 
decreasing in-plane component, which tends to align less 
spins in-plane and thus creates less domain walls, result- 
ing in a larger w [23[ . Based on the T-dependence of Ms 
we estimate that the domains may shrink upon cooling 
down to Tc by less than 10%. 

As alluded to earlier, the samples A and B are nom- 
inally the same and yet they display slightly but no- 
ticeably different superconducting transition tempera- 
ture Tc, coercivity He, coherence length ^(0), and the 
Ginzburg-Landau parameter k. As shown in Fig. 1, 
the width w of the domains for samples A and B is also 
measurably different. These are unavoidable sample-to- 
sample variations for SFB's at these small layer thick- 
nesses, and underscore the importance of performing ex- 
periment with w-tuning on a single SFB sample. 

Figure 2(a) shows the hysteresis loops for sample B 
measured in a SQUID magnetometer. The loops were 
measured in a small external magnetic field, H, between 
-90 Oe and 90 Oe, after setting the domain pattern by 
AC demagnetization at angle 9 followed by cooling the 
sample just below T^. The cycling of the small magnetic 
field has no effect on the domain pattern, and no relax- 
ation of magnetization in the superconducting state has 
been observed. Also included in Fig. 2 (a), with the nar- 
rowest hysteresis loop, is the sample with a saturated F 
layer. The width of the hysteresis loop, AM, increases 
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FIG. 2: (Color online) (a) Hysteresis loops measured with 
SQUID magnetometer after demagnetization of sample B 
with 6 = 0°, 15°, 22.5°, 45°, 67.5°, 75°, 90° and after sat- 
uration (from the outermost to the innermost loop, respec- 
tively). Note that loops for 75° and 90° almost overlap, (b) 
G = AMg/AMsat versus H ior 9 = 0° , 15°, 22.5°, 45°, 67.5°, 
75°, 90° (from top to bottom), (c) G = AMe/AMsat (left 
scale) versus w for H = Q (full squares), H = 22.5 Oe (full 
circles), and ff = 80 Oe (full triangles); and Gu = Ug/Uaat 
for sample C, from [23], (right scale) versus w for H = 1 Oe 
(open squares), H = 80 Oe (open circles) and H = 200 Oe 
(open triangles), (d) Gu = Ue/Usat versus H, extracted from 
transport experiment for sample C [231, for 6 = 0°, 10°, 13°, 
15°, 20°, 45° and 90° (from top to the bottom). 



dramatically as a result of demagnetization indicating a 
large enhancement of vortex pinning by the magnetic do- 
mains. The largest AM is at 6* = 0°, when w is small, 
and decreases as w grows with increasing 9. Even at 9 
= 90° with the largest w of about 0.7 /im, AM is still 
much larger than that of the saturated F layer. We have 
also observed a slight asymmetry in the hysteresis curve, 
particularly at large 9. 

To estimate the enhancement of Jc from the global 
magnetometry results one needs a model of the critical 
state with a specific dependence of Jc on the magnetic 
induction B. The simplest is the Bean model, which 
assumes Jc to be independent of B, leading to the well- 
known prediction of Jc ^ AM |23]. Under this model, 
the ratio AMg/AMgat = G reflects the enhancement of 
Jc induced by the domain pattern as set by the demagne- 
tization, where AMg and AM sat are the hysteresis loop 



width for sample demagnetized at and that of the sat- 
urated F layer respectively. The dependence of G on iJ 
is shown in Fig. 2(b) on a semi-logarithmic scale, omit- 
ting the same dependence for negative H . The G{H) 
dependencies are similar for various 9. The largest value 
of G « 16 is observed for small 9 in the intermediate 
field range of 7? ^ 15 — 30 Oe. The value of G at all 
fields decreases with increasing of 6*, but still retains a 
high peak value of G « 9 at « 90°. Fig. 2(c) shows 
the dependence of G (left scale, full data points) on the 
actual domain width w for three representative field val- 
ues, at the peak {H = 22.5 Oe), low-iJ {H = 0), and 
high-iJ {H = 80 Oe). The G is reduced almost lin- 
early with increasing w at low-_ff (bottom curve). At 
the peak (top curve) the G remains large for narrow do- 
mains, and starts to decrease rapidly for w exceeding 
about 0.45 ^m. The high-_ff behavior (middle curve) is 
intermediate, with the weak suppression of G for small 
w, and more rapid suppression for large w. 

To obtain mo del- independent assessment of Jc, we 
need microscopic measurement of B at various locations. 
For this purpose we employ a linear array of miniature 
Hall sensors to probe locally the dependence of B on the 
distance x from the sample edge. The local magnetic field 
is defined as Hioc = B — fJ-oH, and Jc is obtained from 
the relation ^qJc ~ 2dB/dx [29|. From sample A, we cut 
a 240 iim wide strip and placed a line of 10 sensors, each 
of which 5x5 /um^ in size and 20 /im apart, across the 
strip. An additional sensor, placed a few mm outside the 
sample edge, is used to measure the actual applied field 
as shown in Fig. 3(a). The B{x) is registered simulta- 
neously by Hall sensors while H is swept from up to 
-hlOO Oe, and from -MOO Oe to -100 Oe. The evolution 
of Hioc{x) during the second part of the sweep is shown 
in Figs. 3(b)-(c) for 9 = 0° and 90°. As H is swept, 
the flux remains trapped inside, so that Hioc increases 
in the sample center, eventually reaching a maximum of 
41 Oe (6* = 0°) or 27 Oe {9 = 90°) just after H changes 
sign to negative, while in the saturated sample 5.5 Oe is 
observed. The calculated Jc has a maximum at approxi- 
mately the same H, reaching about 1.5 x 10^ Am^^ for 
9 = 0°, and 8.4 x lO"^ Am^^ for 6* = 90°, while it is only 
1.6 X 10^ Am^^ in the saturated sample. This confirms 
the conclusion of strong pinning enhancement, particu- 
larly for narrow domains at 9 = 0°. Interestingly, Figs. 
3(b)-(c) reveal that the change of Hioc is much smoother 
during the sweep from -1-100 Oe to than during the 
sweep from to -100 Oe, where there are many sudden 
and large drops of Hioc across the whole sample. Such 
abrupt decrease of Hioc indicates the abrupt annihilation 
of flux over the large sample area. It is most likely trig- 
gered by the strong interaction between the positive flux 
still trapped inside and negative flux entering the sam- 
ple. This effect likely contributes to the asymmetry of 
hysteresis seen in Fig. 2(a). 

Figure 3(d) shows the ratio of Jc to J^"* (where J^"* 
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FIG. 3: (Color online) (a) Placement of the sample, with Nb 
adjacent to Hall sensors, (b-c) Hiod^, H) for sample A during 
the field sweep from +100 Oe to -100 Oe for 9 = Q° and 
6 = 90°. The edge and the center of the sample strip are 
at X = 0, and x = 120 /im, respectively, (d) Jcj Jl°'^ vs H 
calculated for a; = 110 /im (x shown by dashed line in (b-c)). 



is Jc in the saturated sample), calculated for x = 110/im, 
the approximate position at which dB j dx is the largest 
(along the dashed line in Figs. 3(b)-(c)). The calcula- 
tion has been made for H > 0, since for negative H Jc is 
highly irregular because of the jumps of Hioc- As shown 
in Fig. 3(d), Jc/Jc""*^ is higher for 6* = 0° in the whole 
i/-range than that for 9 = 90°, the same conclusion as 
that in Fig. 2(b), although the shape of Jc/ Jc""^ , while 
generally similar to G{H) in Fig. 2(b), differs noticeably 
in details. This is not surprising because the local mea- 
surements using an array of Hall sensors are more suscep- 
tible to local irregularities of the flux penetration inside 
the narrow sample, while the global magnetometry mea- 
surement integrates over the large sample. Nevertheless, 
the most essential conclusion that pinning is strongest 
for narrow domains remains unchanged. Similar result is 
obtained at lower T, but the maximal magnitude of the 
Jc/J^* decreases, by about 17% at T/Tc = 0.82, and by 
30% at T/Tc = 0.78. 

The results clearly establish that the strongest pinning 
occurs at the smallest w. This is due to several factors. 
The magnitude of the stray field at the domain center 
scales as 1/w |30j . and moreover, the density of pinning 
centers also increases with decreasing w. As a result, at 
low H the vortices are most effectively pinned when w 
is the smallest. An additional argument comes from the 
comparison of w to the penetration depth, since this is 
the length scale of the vortex-domain interaction. The 
thickness of Nb in sample B is i = 76 nm, slightly larger 
than A(0) « k£,{0) ~ 64 nm. As T increases, A(r) grows 
as - <(0)/\/l - T/Tc, so at T/Tc = 0.91 we obtain 
A « 220 nm. However, at high T one approaches the 
limit of A ^ t, for which one should use the effective 



penetration depth, A = A^/i w 640 nm [27|. While this 
higher hmit is not quite reached here, it is hkely that A 
approaches, and perhaps even exceeds, the average do- 
main width of w w 310 nm for 6 — 0°. Since in the 
quasi-ordered domain pattern there are areas in which 
w is smaller than average, vortices cannot move freely 
along domains, thus enhancing the pinning. The value 
of w grows with increasing 9, and the vortices flow more 
easily, thus diminishing the value of G. 

As shown in Fig. 2(b), for each value of 9, G increases 
with H reaching a maximum presumably when the vortex 
lattice becomes commensurate with the magnetic domain 
pattern. This is followed by a decrease of G, when the 
vortex density exceeds that of the domain-induced pin- 
ning centers. Since the domain density decreases with 
increasing w, the suppression of G at high H is weaker 
for smaller w as shown by the blue triangles in Fig. 2(c). 

The commensurability of vortex pinning can be ad- 
vantageously observed using an ordered array of pinning 
centers, which lead to enhanced pinning at well-defined 
matching fields \M,- In the SFBs the vortices are con- 
fined to the irregular domains of one sign. In the case of 
regular stripe domain pattern with the period 2w single 
or multiple chains of vortices are pinned by each domain, 
as recently shown by scanning tunneling microscopy 22 1 . 
We have recently performed transport measurements on 
sample C demagnetized at different 9 to capture differ- 
ent domain widths and we have inferred a similar chain- 
pinning phenomenon [23| . At T sufficiently below Tc the 
resistance is thermally activated, with the activation en- 
ergy for vortex pinning, U, enhanced by magnetic do- 
mains. The U{H) exhibits maxima at matching fields 
when narrow domains (small 9) pin single vortex chains, 
and wider domains (large 9) pin double vortex chains. 

It is interesting to compare the activation energy U 
from transport measurements on sample C [23] with the 
present data on sample B. From the transport results, we 
calculate the enhancement of U induced by demagnetiza- 
tion at T/Tc = 0.91, using the parameter Gu — Ug/Usat, 
where Ug is for sample C demagnetized at 0, and Usat 
is for saturated sample. The values of Gu versus H for 
several values of 9 are shown in Fig. 2(d). One notes 
that Gu and G display strikingly similar H dependence. 
Furthermore, both display the largest value at the small- 
est values of 9, decreasing monotonically with increasing 
6. These systematic dependencies are illustrated in Fig. 
2(c), where G (filled symbols) and Gu (open symbols) are 
for three representative fields, from the low field, peak, 
and high field regions {H = 0, 22.5, and 80 Oe for G, 
and i? = 1, 80, and 200 Oe for Gu, respectively). The 
qualitative similarity suggests a common origin, which is 
pinning of vortex chains by domains. 

However, there are also noticeable differences between 
G{H) and Gu{H). One notes that the peak in Gu{H) 
has a narrower width but it is located at a higher H 
than that of G{H), whereas G reaches a maximum mag- 



nitude about 5 times larger than that of Gu- Most of 
these differences are due to the different Nb thickness in 
sample B (76 nm) for G{H) and sample C (20 nm) for 
Gu (H) . In addition, sample B shows larger dispersion of 
w. The dispersion of lu contributes directly to the broad- 
ening of matching field, while a thicker Nb layer leads 
to a progressive reduction of the domain-induced stray 
field (by up to about 25% on the other side of Nb [SOj), 
and therefore reduction of the density of pinning centers 
across t, shifting G{H) peak to smaller H. The difference 
in magnitude between G and Gu is also mainly related 
to thickness. The strength of the magnetic interaction 
roughly depends on the ratio of the domain width w, to 
the range of magnetic interaction, given approximately 
by A. In the thicker sample B the magnetic interaction 
is stronger because w is comparable to A, whereas in the 
thinner sample C, A is about 4 times larger so that the 
magnetic interaction is substantially reduced. 

In conclusion, we have demonstrated a method to in- 
duce large enhancement of vortex pinning in the su- 
perconductor/ferromagnet bilayers using a ferromagnetic 
layer with perpendicular magnetic anisotropy. By de- 
magnetizing the sample at an angle to the sample surface, 
we obtain domain patterns with equal +/— domains but 
different domain width w. This unique attribute allows 
a single bilayer to acquire different domain width to ex- 
ert tunable vortex pinning on a superconducting layer. 
Magnetometry and local measurements using an arrays 
of Hall sensors show conclusively vortex pinning enhance- 
ment, by a factor of more than 10 at domain width w of 
about 310 nm, much larger than previously obtained. 

This work was supported by Polish MNiSW grant 
N202 058 32/1202, by NSF grant DMR05-20491, by the 
French-Pohsh Bilateral Program PICS 4916, and by the 
European Union within the European Regional Devel- 
opment Fund, through the Innovative Economy grant 
POIG.01.01.02-00-108/09. 



Current address: Materials Science Division, Argonne 

National Laboratory 

M. Velez, J.I. Martn, J.E. Villegas, A. Hoffmann, E.M. 

Gonzlez, J.L. Vicent and Ivan K. SchuUer, J. Mag. Magn. 

Mat. 320, 2547 (2008). 

I. F. Lyuksyutov and V. L. Pokrovsky, Adv. Phys. 54, 

67 (2004). 

A. Yu. Aladyshkin, A. V. Silhanek, W. Gillijns, and V. 

V. Moshchalkov, Supercond. Sci. Technol. 22, 053001 

(2009). 

L. N. Bulaevskii, E. M. Chudnovsky, and M. P. Maley, 

Appl. Phys. Lett. 76, 2594 (2000). 

Yu. I. Bezpyatykh, W. Wasilevski, M. Gajdek, I. P. 

Nikitin and S. A. Nikitov, Sov. Phys. -Solid State 43, 

1827 (2001). 

M. V. Milosevic, S. V. Yampolskii, and F. M. Peeters, 

Phys. Rev. B 66, 174519 (2002). 



[7] S. Erdin, I. F. Lyuksyutov, V. L. Pokrovsky, and V. M. 

Vinokur, Phys. Rev. Lett. 88, 017001 (2002). 
[8] E. B. Sonin, Phys. Rev. B 66, 136501 (2002). 
[9] R. Laiho, E. Lahderanta, E. B. Sonin, and K. B. Traito, 

Phys. Rev. B 67, 144522 (2003). 
[10] M. A. Kayah and V. L. Pokrovsky, Phys. Rev. B 69, 

132501 (2004). 
[11] S. Erdin, Phys. Rev. B 73, 224506 (2006). 
[12] A. Garcia-Santiago, F. Sanchez, M. Varela, and J. Te- 

jada, Appl. Phys. Lett. 77, 2900 (2000). 
[13] X. X. Zhang, G. H. Wen, R. K. Zheng, G. C. Xiong, and 

G. J. Lian, Europhys. Lett. 56, 119 (2001). 
[14] D. B. Jan, J. Y. Couher, M. E. Hawley, L. N. Bulaevskii, 

M. P. Maley, Q. X. Jia, B. B. Maranville, F. HeUman, 

and X. Q. Pan, Appl. Phys. Lett. 82, 778 (2003). 
[15] M. Lange, M. J. Van Bael, V. V. Moshchalkov, and Y. 

Bruynseraede, AppL Phys. Lett. 81, 322 (2002). 
[16] M. Z. Cieplak, X. M. Cheng, C. L. Chien, and Hai Sang, 

J. Appl. Phys. 97, 026105 (2005). 
[17] M. Z. Cieplak Z. Adamus, A. Abal'oshev, I. Abal'osheva, 

M. Berkowski, X. M. Cheng, Hai Sang, C. L. Chien, Phys. 

Status Solidi C 2, 1650 (2005). 
[18] M. Feigenson, L. Klein, M. Karpovski, J. W. Reiner, and 

M. R. Beasley, J. Appl. Phys. 97, 10J120 (2005). 
[19] V. Vlasko-Vlasov, U. Welp, G. Karapetrov, V. Novosad, 

D. Rosenmann, M. lavarone, A. Belkin, and W.-K. Kwok, 

Phys. Rev. B 77, 134518 (2008). 



[20] A. Belkin, V. Novosad, M. lavarone, J. Pearson, G. Kara- 
petrov, Phys. Rev. B 77, 180506 (2008). 
[21] V. K. Vlasko-Vlasov, U. Welp, A. Imre, D. Rosenmann, 

J. Pearson, and W. K. Kwok, Phys. Rev. B 78, 214511 

(2008). 
[22] G. Karapetrov, M. V. Milosevic, M. lavarone, J. Fedor, 

A. Belkin, V. Novosad, and F. M. Peelers, Phys. Rev. B 

80, 180506(R) (2009). 
[23] L. Y. Zhu, M. Z. Cieplak, and C. L. Chien, Phys. Rev. B 

82, 060503(R) (2010). 
[24] V. Vlasko-Vlasov, U. Welp, W. Kwok, D. Rosenmann, 

H. Glaus, A. A. Buzdin, and A. Melnikov, Phys. Rev. B 

82, 100502(R) (2010). 
[25] M. lavarone, A. Scarfato, F. Bobba, M. Longobardi, G. 

Karapetrov, V. Novosad, V. Yefremenko, F. Giubileo, 

and A. M. Cucolo, Phys. Rev. B 84, 024506 (2011). 
[26] C. Visani, P. J. Metaxas, A. CoUaudin, B. Calvet, R. 

Bernard, J. Briatico, C. Deranlot, K. Bouzehouane, J. E. 

Villegas, arXiv:1107.1122i^l. 
[27] M. Tinkham, Introduction to superconductivity, Dover 

Publications, 2004. 
[28] C. P. Bean, Phys. Rev. Lett. 8, 250 (1962); Rev. Mod. 

Phys. 36, 31 (1964). 
[29] E. H. Brandt, Phys. Rev. B 54, 4246 (1996). 
[30] G. M. Maksimova, R. M. Ainbinder, D. Y. Vodolazov, 

Phys. Rev. B 78, 224505 (2008). 



